Solitonic Excitations in Linearly Coherent Channels of Bilayer Quantum Hall Stripes 
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In some range of interlayer distances, the ground state of the two-dimensional electron gas at filling 
factor v = AN + 1 with N = 0,1, 2, ... is a coherent stripe phase in the Hartree-Fock approximation. 
This phase has one-dimensional coherent channels that support charged excitations in the form of 
pseudospin solitons. In this work, we compute the transport gap of the coherent striped phase due 
to the creation of soliton-antisoliton pairs using a supercell microscopic unrestricted Hartree-Fock 
approach. We study this gap as a function of interlayer distance and tunneling amplitude. Our 
calculations confirm that the soliton-antisoliton excitation energy is lower than the corresponding 
Hartree-Fock electron-hole pair energy. We compare our results with estimates of the transport gap 
obtained from a field-theoretic model valid in the limit of slowly varying pseudospin textures. 

PACS numbers: 73.43.-f, 73.21. Fg, 73.20.Qt 



I. INTRODUCTION 

It is well known that the ground state of the two- 
dimensional electron gas (2DEG) in single layer quantum 
Hall systems near half-odd integer filling factors in Lan- 
dau levels N > 2 i.e. for v = 9/2, 11/2, ... is a striped 
state responsible for a strong anisotropy in the conduc- 
tivity tensor of the 2DEG. This state was predicted on 
the basis of Hartree-Fock calculations^ and has been ex- 
tensively studied experimentally^ 

When the interlayer distance, d, in a bilayer quantum 
Hall system at filling factor v is large, one expects the sys- 
tem to behave as two isolated two-dimensional electron 
gases (2DEG) with filling factor v/2. It is then natural 
to infer that the ground state of the 2DEG in a bilayer 
should be a striped state at v = 4N + 1 at sufficiently 
large interlayer distances. On the other hand, it is known 
that, at v = 4N + 1 interlayer interactions can lead to 
a homogeneous ground state with spontaneous phase co- 
herence between the layers when the interlayer distance 
is comparable with the separation between electrons in 
a single layer. One might then conjecture that, as the 
interlayer separation is decreased, the striped state ac- 
quires a certain degree of coherence due to the interlayer 
interaction. This conjecture was first studied by Brey 
and Fertig^ who showed that, as d is increased from zero 
the bilayer ground state goes from a uniform coherent 
state (UCS) at small interlayer separations to a coherent 
striped phase (CSP) at d > d± and then into a modulated 
striped state (or anisotropic Wigner crystal) at d > di. 
The interlayer coherence is lost in the modulated stripe 
state. The range [di,da] increases with 

The coherent striped phase shown in Fig. ^is a state 
where charge density waves in the two layers are shifted 
by £/2 where £ is the period of the stripes in one layer. 
The most interesting aspect of the CSP is that in the 
regions where the charge densities in both layers "over- 
lap" (in the plane of the two-dimensional electron gas 
(2DEG)), the electrons are effectively in a linear superpo- 
sition of states of the form = (\R) + \L)) /V2 where 



R, L indicates the right and left wells. The interlayer co- 
herence is then maintained but only along linearly coher- 
ent regions (LCR's) whose width decreases as d increases. 
The CSP is most easily represented in the pseudospin 
language where an up (down) pseudospin is associated 
with the right (left) well. The CSP is a pseudospin den- 
sity wave where the pseudospins oscillates in the xz plane 
and the LCR's are the one-dimensional regions where the 
pseudospins lie along the x direction in the xy plane. 




FIG. 1: Guiding center density in the right (dark surface) and 
left (light surface) wells and pseudospin pattern in the coher- 
ent stripe phase. The arrow indicates one linearly coherent 
channel (LCR). 

In a previous work^, we have computed the collective 
excitations of the CSP and showed that the low-energy 
modes of this phase could be described by an effective 
pseudospin wave hamiltonian. We have also shown 6 that 
the application of a parallel magnetic field gives rise 
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to a very rich phase diagram for the 2DEG involving 
commensurate-incommensurate transitions with distinc- 
tive signatures in the collective excitations and tunneling 
I — V. A very exhaustive study of the phase diagram of 
the 2DEG in the presence of a parallel magnetic field, in 
higher Landau levels, has also been published by Daw- 
Wei Wang et al£&. 

The band structure of the CSP is shown in Fig. [21 
In the Hartree-Fock approximation, the energy gap of 
this system corresponds to the excitation of an electron- 
hole pair in a coherent channel (a pseudospin flip in 
the xy plane) and is finite if the tunneling parameter 
( ^ 0. An estimate of this gap, taking into account 
some quantum fluctuations, has been done by E. Papa 
et al&. However, Brey and Fertigi pointed out that, in 
analogy with spin (pseudospin) skyrmion excitation in 
single (double) layer quantum Hall systems at v = 1, 
the lowest-energy charged excitation should be a pseu- 
dospin soliton (or antisoliton) in a coherent channel and 
the gap should be given by the energy required to cre- 
ate a soliton-antisoliton pair. A pseudospin soliton of 
charge q — e corresponds to a 2tt rotation of the pseu- 
dospin in the xy plane. As for skyrmions or bimerons, 
the size of these solitons is determined by a competition 
between tunneling energy (which favors small solitons) 
and interwell exchange energy and Coulomb interaction 
which favors slowly varying pseudospin textures (large 
solitons). 

In this work, we compute the energy gap of the CSP 
due to the excitation of a soliton-antisoliton pair as a 
function of tunneling and interlayer distance. We use 
a supercell microscopic unrestricted Hartree-Fock ap- 
proach to extract the energy of a single soliton from that 
of a crystal of solitons localized in the LCR's at filling 
factor v = 47V + 1 + Av. Our calculation shows that 
a soliton-antisoliton pair has a lower energy than the 
electron-hole pair so that these topological excitations 
will be important in determining the transport properties 
of the CSP. For completeness, we also compute the energy 
gap of the CSP using a simple field-theoretic model based 
on the sine-Gordon Hamiltonian where an exact solution 
for the pseudospin soliton can be obtained. This model 
does not contain all the terms included in the microscopic 
approach, but, for slowly varying pseudospin textures, it 
should give a fair estimate of the energy gap. We actu- 
ally improve on this model by taking into account that 
the channels have a width that depends on the interlayer 
distance d and also by taking into account the interaction 
of the pseudospins in different channels and the Coulomb 
interaction between different portions of the topological 
charge densities. 

The paper is organized as follows. In Sec. II, we de- 
scribe the phase diagram of the 2DEG in the bilayer sys- 
tem at filling factors v = 47V +1 and v = 47V + 1 + Az; and 
define the domain of existence of the soliton crystal from 
which we want to compute the soliton energy. In Sec. 
Ill, we introduce the simple field-theoretic model and the 
exact solution for the pseudospin sine-Gordon solution. 



Section IV discusses the supercell method that we use to 
extract the energy of a single soliton from that of a crys- 
tal of solitons. The removal of the soliton-soliton energy 
is discussed in Sec. V. Section VI discusses our numerical 
results. We conclude in Sec. VII. Details of the deriva- 
tion of the microscopic expression for the parameters of 
the field-theoretic model arc given in the appendix. 
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FIG. 2: Band structure of the coherent stripe phase. The 
greyed states represent filled states at v — 47V + 1. The 
Hartree-Fock gap is also indicated. It corresponds to the ex- 
citation of an electron-hole pair in one of the linearly coherent 
channels. 



II. PHASE DIAGRAM OF THE 2DEG AROUND 

v = 47V + 1 

In this section, we review the phase diagram of the 
2DEG at filling factor v — 47V + 1 where the coherent 
striped state is found and at filling factors slightly above 
v = 47V + 1 in order to find the range of interlayer dis- 
tances where a crystal of solitons localized in the LCR's is 
stable. We need the energy of this soliton lattice in order 
to compute the gap energy as we explained in the intro- 
duction. To establish the phase diagram, we compute the 
energy of different electronic phases in the Hartree-Fock 
approximation in order to find the one that minimizes 
the total energy at a given value of v, d, and t. The order 
parameters for the different phases are the expectation 
values of the density operator projected onto the Landau 
level TV of the partially filled Landau level (the guiding 
center density), i.e., 

(P%W) = ^ E e- W+X ' )/2 5x,x>-^ (1) 
v ^ x,x> 

X \ c X,i,N C X',i,x) , 

where i,j are layer indices and X, X' are guiding cen- 
ter coordinates^. We make the usual approximation of 
assuming that the filled levels are inert. We also neglect 
Landau level mixing and assume that the electron gas in 
the partially filled level is fully spin polarized. In a crys- 
tal phase, (/9^(q)) is non zero only for q = G where G 
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is a reciprocal lattice vector of the crystal. Defining the Hartree and Fock interactions 
I 



and 
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where LfJ (x) is a generalized Laguerre polynomial, 
Jo (x) is the zeroth-order Bessel function of the first kind 
and the form factor 



A; 




if i =3, 
if * ^J, 



(4) 



the Hartree-Fock energy per electron at total filling factor 
v = 4N + v can be written as 



E 



(5) 



with 



7 Re [(p B N L (0))] (6) 
h E E H iJ ( N > N > G ) (pn (-G)) 6# (G) 



E E *m W ^> G ) (pn (-G)) (a# (G) 
~f E E ^(n^'.O) 
-- V X iA (n,N,0)u. 

n<N 

In this last equation, N e is the total number of electrons 
in the 2DEG, t is the tunneling strength (in units of 
(e 2 /Ki), with k the dielectric constant of the host mate- 
rial and I = ^Jhc/eB the magnetic length). 

The last two terms in Eq. @) give the interaction be- 
tween electrons in the filled levels and between electrons 
in the filled levels and electrons in the partially filled level 
N. As we will show later, the filled levels contribute to 
the quasiparticle energies, but not to the charge gap. 

The set of (p^(G))'s corresponding to one particular 
electronic phase is found by solving the equation of mo- 
tion for the one-particle Green's function in the Hartree- 
Fock approximation. The method is described in detail 
in Ref. 



The band structure of the CSP contains two bands 
E± (X), as shown in Fig. At exactly v = 4N + 1, 
the lowest-energy band is completely filled and the sys- 
tem is gapped even in the absence of tunneling. In 
fact, in the uniform coherent state that occurs for val- 
ues of d for which stripe ordering had not set in, the 
band structure consists of two straight lines separated 
by a gap Ayes = (2t + 2X R x(N, N; 0)) (e 2 / ' nl) with 
Aucs — ► 2t as d — > oo. In the CSP, the energy bands are 
periodically modulated in space with the maxima (min- 
ima) of the valence (conduction) band at the locations 
of the LCR's. At the Hartree-Fock level, the energy gap 
is the energy needed to excite an electron from a maxi- 
mum of the valence band to a minimum of the conduction 
band. This excitation corresponds to a single spin flip lo- 
calized in one LCR. The decrease in the HF gap in the 
CSP is due not so much to the reduction of Xr^{N, N; 0) 
with d as to the increase in intralayer correlations that 
increases the with of the modulations in E± (X). As d 
increases, the charge modulations get sharper up to the 
point where the stripes become square waves at very large 
d. Correspondingly, the width of the LCR's decreases 
with d since interwell coherence and charge modulation 
compete with each other. 

In analogy with the excitations of skyrmions in single 
quantum well and bimerons in bilayer systems at v = 1, 
Brey and FertigS noted that a lower-energy excitation 
could be achieved by exciting a pseudospin soliton in the 
LCR instead of a simple electron-hole pair. The pseu- 
dospin soliton corresponds to a 2-7r rotation of the pseu- 
dospin in one LCR. A slowly varying pseudospin config- 
uration like that in a soliton has lower exchange energy 
than a single pseudospin flip but the cost in tunneling 
energy is increased. As for skyrmions or bimerons, an 
optimal size for the soliton is obtained at given values 
of v, d and t. The energy cost for this optimal soliton 
should be compared with the Hartree-Fock electron-hole 
pair excitation to determine whether or not these topo- 
logical excitations are energetically favorable. 

In a quantum Hall system, the relation between the 
charge density of the solitons and their pseudospin tex- 
ture (at v — 1) is given by the Pontryagian densityii 
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6 {P (r)> = 8^ £abcSa (r) £ ' ASb (r) djSc (r) ' (?) 

where £y and e a bc are antisymmetric tensors and S (r) 
is a classical field with unit modulus representing the 
pseudospins and S (p (r)) is the guiding-center density. If 
we write a general solution as 

S x (r) = sin0(r)cos(^(r), (8) 
S^fr) = sin 6» (r) sin V3 (r) , (9) 
S z (r) = cos#(r), (10) 

then the induced density takes the simple form 

Sp (r) = sin0 (r) [Vy> (r) x Vtf (r)] • 2. (11) 

In a LCR, the polar angle of the pseudospins 6 = ir/2. 
If a soliton is present in this LCR, then <p (r) rotates 
by ±27r along the channel (oriented in the y direction). 
As discussed below, this is a generalization of a soliton 
in the sine-Gordon modeliS. We also have that, in the 
CSP, V0 (r) 7^ in the LCR's and so the solitons carry 
a charge by virtue of Eq. lfTT|l . 

In the case where pseudospin solitons are the lowest- 
energy excitations of the CSP, we expect that the ground 
state at v = AN + v will be a crystal of solitons localized 
in the LCR's. Table I shows that the range of interlayer 
distances where the CSP is the system's ground state at 
v = 4N+ 1 increases with the Landau level index. In this 
work, we choose to study the phase diagram in Landau 
level N = 2. We show in Fig. EJthe energy per electron 
for different electronic phases in N = 2 as a function of 
interlayer distances and for three values of the tunneling 
parameter t = 0,0.01 and 0.06. The filling factor is v = 
9.2. The contribution from the filled levels is not included 
in this calculation since it depends only on v and is thus 
the same for all phases. At small interlayer distances, 
where the ground state at v — 9 is a UCS, the ground 
state at v = 9.2 is a one-component hexagonal Wigner 
crystal (HWC). In this phase, a crystal of electrons of 
pseudospin S x = —1/2 and filling v = 0.2 sits on top of 
a liquid of pseudospins S x = +1/2 and filling 9.0. There 
is no pseudospin texture in that state and, in particular, 
no bimerons in contrast with the situation in the lowest 
Landau leveti^ where the ground state is a crystal of 
bimerons. In fact, we find that bimeron excitations are 
not relevant in N — 2 even in the limit of vanishing t. 

For interlayer distances where the CSP is found at 
v = 9, the ground state of the 2DEG at v = 9.2 is a 
centered crystal of pseudospin solitons localized in the 
LCR's. We note that there are many possible choices for 
the lattice structure of this crystal, since solitons may 
or may not be present in every LCR, depending on the 
commensuration of the lattice of solitons and the under- 
lying stripe state, and it is likely that there are phase 



Landau level <h/l <hjj. Djl 






1.2 


1.65 


0.45 


1 


0.8 


1.45 


0.65 


2 


0.6 


1.6 


1.00 



TABLE I: Critical interlayer distances d\jl and d^/l at t = 
for the transition UCS-CSP and CSP-modulated striped 
state. The last column gives the range of interlayer distances 
D/£ = dill — di/i for which the CSP is the ground state in 
Landau level N. 



transitions among these different states as the filling fac- 
tor is varied. For the choice of parameters in this study, 
the lowest energy state has solitons in every channel. We 
found however that a similar state with solitons in ev- 
ery second channel but with the same filling factor has 
very nearly the same energy. Figure 0] shows an example 
of the charge distribution as well as the pseudospin tex- 
ture associated with a centered rectangular soliton crys- 
tal. Since the focus of this study is on the energetics of 
single solitons, we will use only the structure illustrated 
in Fig. 0]for our quantitative analysis below. 

At large interlayer distances, we find that the ground 
state of the 2DEG at v = 9.2 is a superposition of 
two shifted triangular bubble crystalsi with partial fill- 
ing v = 0.6 in each well. Because v > 0.5 the bubbles 
are clusters of holes and not electrons. We find that the 
number of holes per bubble is M = 3 in agreement with 
previous Hartree-Fock calculation in single quantum well 
systems^. 

III. FIELD-THEORETIC MODEL 

We use two different approaches to compute the energy 
gap due to the excitation of soliton-antisoliton pairs. The 
first one is a field-theoretic calculation valid in the limit 
of slowly varying pseudospin textures. It is explained in 
this section. The second one is a microscopic approach 
where the energy of one soliton is computed from that 
of a crystal of solitons by removing the soliton-soliton 
interaction. We call this method the supercell approach. 
In principle, this second method is not restricted to small 
gradient of the pseudospin texture and includes terms 
neglected in the field-theoretic model. We expect it to 
be more accurate than the field-theoretic approach. 

In the field-theoretic approach, we evaluate the en- 
ergy to create a pseudospin soliton by making a long- 
wavelength expansion of certain terms in the Hartree- 
Fock Hamiltonian. We follow the procedure developped 
in details in Ref. To keep the discussion as brief as 
possible, we give here only the main results of this model. 
Full details are provided in the appendix. 

There are three main contributions to the energy 
needed to create a pseudospin texture in a LCR. Since in 
the ground state the in-plane pseudospin component in 
a LCR is fully polarized along S x , adding a pseudospin 
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FIG. 3: Hartree-Fock ground state energy per electron as a 
function of interlayer distances at filling factor v — 9.2 and 
for (a) t = 0; (b) t = 0.01; (c) t = 0.06. The vertical lines 
indicate the position of the phase transitions. 
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FIG. 4: Representation of the soliton crystal at d/l = 1.2, £ = 
0.01 and i> = 9.1. The distance between two solitons in a chan- 
nel is a. (a) Guiding-center densities pun (x, y) , pll (x, y) and 
p (x, y) = pfli? (x, y) + pll (x, 2/) at y = 0; (b) pseudospin tex- 
ture showing the solitons localized in the channels. 



texture has a tunnel energy cost when t ^ because of 
the interaction of the texture with the other channels. A 
second contribution comes from the interlayer exchange 
interaction which is responsible for the pseudospin stiff- 
ness p s . As we mentioned above, the exchange interac- 
tion favors pseudospin textures that vary slowly in space. 
A third contribution must be considered in our model in 
order to get agreement with the microscopic approach. 
It is the Coulomb interaction between different portions 
of the soliton in a channel. This interaction favors large 
solitons. 

If the coherent channels are oriented along y and are 
considered as effectively one-dimensional, then the en- 
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ergy cost to make a pseudospin texture on top of the 
ground state where all pseudospins point in the x direc- 
tion in each channel is 



SE = dy 



V dy J 



(12) 



where (p(y) in the azimuthal angle of the pseudospins. 
Eq. H12fl is valid if we ignore the third contribution men- 
tionned above. The parameters p s and T are the effective 
stiffness and tunneling parameters. These parameters de- 
pend on the precise shape of the LCR's as well as on the 
interaction between pseudospins of different channels. In 
the appendix, we derive a microscopic expression for each 
of these parameters in terms of the order parameters of 
the CSP. We show that the effective stiffness is given by 



-1 



16tt 2 £ 2 



^) J dq x \tt(q x ) 



d 2 X RtL (N,N;q) 



dql 



(13) 



where 



Or 



sin [(Gs - 
(G x -q x )£/4 



(14) 



is a form factor that takes into account the shape of the 
channel centered at x — 0. Also, £ is the interstripe dis- 
tance in the CSP, G x = 2irn/£ with n = 0, ±1, ±2, ... and 

(p%(G x )) — Re /p^r' L (G x )^ . If we define the parame- 
ter G x = 4irn/£ and 



Jx(q) = -X fl>£ (N,N;(0, 



then the parameter T can be written as 



(15) 



T = 



2tt£ 2 



m(q x = o) - j x (g„o) \n(G x )\ 2 



(16) 



The second and third terms in Eq. (|16fl come from the 
fact that, because of the pseudospin stiffness, there is 
an energy cost to rotate the pseudospins in one channel 
when the pseudospins in the other channels remain fixed 
in their ground state position. The contribution of these 
two terms increases the effective tunneling strength T. 
Since the energy cost to create a pseudospin soliton is 
given by E s = 8Vp s T we see that this second term keeps 
E s finite even when t = 0. 

In this field-theoretic model, the energy to create an 
antisoliton is the same as that needed to create a soliton 
and the charge gap is simply given by 

A = 16y/p/T. (17) 

From the energy functional of Eq. (|12|) , we get that the 
static solution that minimizes the energy must satisfy the 
sine- Gordon equation 



— — = — smip{y). 



(18) 



The sine-Gordon (or kink) soliton is a solution of this 
equation. It is given by 



<p(y) = 4 tan 



(19) 



The length of the soliton can be defined as 



(20) 



With the energy functional of Eq. I|12|) , we find numer- 
ically that both p s and T decrease rapidly with d but the 
size of the soliton L s decreases with increasing d. This 
behavior is opposite to what we obtain in the microscopic 
calculation where the soliton size increases with d. As we 
mentionned above, it is necessary to include the Coulomb 
interaction between different part of the solitons in order 
to get the soliton size to increase with d. This leads to 
the term (full details are given in the appendix) 



5E, 



Coul 



32tt 2 



dy 



dy — Yeff{y 



y') 



dy' 



. (21) 

Inclusion of this term in in the energy functional intro- 
duces a nonlocal non-linear term in the differential equa- 
tion for the soliton and the resulting equation is very dif- 
ficult to solve. Following S. Ghosh and R. Rajaramaolfii 
who use a similar procedure in their calculation of the 
energy of CP 3 skyrmions in bilayers, we make the fol- 
lowing approximation. We insert a pseudospin texture 
9?(y) = 4 tan -1 [e _J, / L =] into the total energy functional 
including the Coulomb integral and evaluate is as a func- 
tion of L*. We then minimize the total energy with re- 
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spect to the length L* to obtain the energy and length 
of the soliton. In this way, we find a soliton length that 
increases with d as in the microscopic approach. The 
procedure is described in details in the appendix. 



IV. THE SUPERCELL MICROSCOPIC 
HARTREE-FOCK METHOD 

Let £csp be the energy per electron in the CSP at 
v = AN + 1 and magnetic field Bo in units of e 2 / 'k£q. If 
the number of electrons is kept constant and the magnetic 
field is decreased (to B{) or increased (to B 2 ) such that 
the filling factor becomes v = AN±v, then a finite density 
n qp — \v — 1| /2tt£1 2 of quasiparticles (solitons for v > 1 
and antisolitons for v < 1) are created in the CSP. At zero 
temperature, we expect these quasiparticles to crystallize 
and to be localized in the LCR's of the CSP. In the limit 
where only one quasiparticle is created (y — ► 1), we can 
define the quasiparticle energy as 



lim • 



£sc 



--CSP 



k£ 



(22) 

where esc is the energy per electron in the soliton crystal 
(SC) in units of e 2 / 'k£ with N qp solitons and E+ p (E~) 
is the energy to create one soliton (antisoliton). 

The quasiparticle energy defined in this way, with the 
number of electrons kept constant, is refered to as the 
"proper" quasiparticle energy by Morf and Halperm 16 . 
Other definitions are also possible. For example, the 
"gross" quasiparticle energies (or chemical potentials) are 
defined by 



E{N e = N^ + l)-E{N e 



H~ = E(N e = N^) — E (N e = N^ — l), 



(23) 
(24) 



where N$ is the degeneracy of the Landau levels at a 
magnetic field Bq such that v = AN + 1. The energy 
E {Ncj,) is the total energy of the CSP, and E (N</, ± 1) 
is the total energy of the CSP with one more (less) par- 
ticle in the form of a soliton (antisoliton). In this case, 
the magnetic field is kept constant while the number of 
particles changes. At zero temperature, this is precisely 
the definition of the chemical potentials at filling factors 
slightly above or below v — AN + 1. 

The different definitions of the the quasiparticle ener- 
gies lead to different numerical values. As discussed by 
MacDonald and Girvinii, however, the numerical value 
of the gap, A, is the same for both definitions so that we 
can write 



A 



/' 



e: 



(25) 



With the formalism described in Sec. II, we can easily 
compute the Hartree-Fock energy of a crystal of solitons 
located in the coherent channels of the bilayer. That is, 
we can compute £sc: nn d E^ p and then the energy gap. 
However, there are several difficulties with this method 



that we address in this paper. The first one is that the 



limit 



1 cannot be achieved numerically since that 



would require infinite matrices in the equation of motion 
for the single-particle Green's function. In this work, we 
have succeeded in computing esc at filling as small as 
v = 1 ± 0.02. The second difficulty is that, when a fi- 
nite density of quasiparticles is present, esc includes the 
interaction energy between quasiparticles. This interac- 
tion energy must be computed and removed from esc- 
A third difficulty is related to the size of the solitons 
(antisolitons). In Sec. Ill, we saw that the soliton size 
becomes very large when the tunneling energy t — > or 
when d is large. In this case the size of the soliton is not 
given by Eq. H20fl but is limited by the lattice constant of 
the soliton crystal. The quasiparticle energy, then, can- 
not be computed reliably when the tunneling term is too 
small or the interlayer distance too big. 

We now describe in more details our evaluation of E^p. 
To avoid computing numerically the energy of the anti- 
soliton crystal as well as that of the soliton crystal, we use 
the particle-hole symmetry of the Hamiltonian around 
v = AN + I to relate the energies of the two crystals with 
the same filling of quasiparticles. We define state as 
the CSP at v — AN + I, state 1 as the soliton crystal at 
v\ = AN + v\ and state 2 as the crystal of antisolitons at 
V2 = AN + V2- The filling factors v 2 = 2 — v\ so that the 
lattice constants a\ and a 2 of the two crystals are related 
by £\/ai — £2/0-2- The Hartree-Fock energy per electron 
of the three states are given by Eq. JSJ which we rewrite 
here as 



E m 

N K 



— J e m {v m ) + —A (v m ) 



We have defined 



2t r 

Em (^rn,) = -—Re (/9^' L (0)) T 



(27) 



Zj is <> y 1 ... • * TTi » 



i,j G^O 



^EE X « (pn (&' ( G ) 

Zit/ity-i \ I 771. \ 



id G 

which is the energy per electron in the partially filled 
level. The last term in Eq. (|26|1 is the interaction energy 
with the filled level with 



K(vi) = -2Ai - A 2 ^i, 



where 



Ai = E X R , R {n,n',0), 

n<N n'<N 

A 2 = X i;i (n,N,0). 



(28) 

(29) 
(30) 



n<N 



From Eqs. (|23|l and l|24|) . it is easy to see that the 
cyclotron and Zeeman energies do not contribute to the 
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transport gap A and so can be ignored in Eq. (|26|l . This 
is also true of the filled levels since their contribution to 
the quasiparticle energies are given by 



(E + ) 
v qpl f.i. 



lim i - ,i , 

JVgp-fi \v\ — 1| \nt\ ) v\ 



— lim 
N qp ^i \vi - 1 AN 



A fa) (31) 

2 

k£ , 



' I A(1) / ' 



e 



1 



-A 2 + 3Ai 



and 



E 



ip) f.i. 



lim — — - — r 

N gp ^l \v 2 - 1 



lim 



1 

V 2 



A fa) 
1 



(32) 



/ N qp ^i \v 2 - 1| 4AT + 1 







(«4>) 


^A 2 + 3A! 





so that ( E2 



"qp) f.i. V^qp) f.i 
equations, we have used 



e 

k£ 



A(l) 

0. In deriving these two 
X (33) 



From the electron-hole symmetry, we get 



£2 



El 



?1 



X(0) 



where 



A^(0) = X fl)fl (iV,7V,0) 



(34) 



(35) 



Note that Eq. I)34|l is exact only in the limit where N qp 
1 because the inter-well Hartree and Fock interactions 
contained in e m depend on the ratio d/i and we have 

d/ti £ d/e 2 . 

Combining all results, we have for the energy gap 



A = lim — — v\ 

Ai/->0 Av 

1 

+ lim — - 

Ai/^0 Av 



where we have defined 

£csp = 



X (0) - 2e CS p 



e 

k£ 



e 

k£ 



1 



AN+1 



£csp- 



Simplifying, we get finally 



A = lim 

Ai/^0 



2-^-ei -^-scsp + X(0) 
Av Av 



(36) 
,(37) 

(38) 



TH- (39) 



e 

k£ 



We remark that the change in the magnetic length £ due 
to the change in the magnetic field makes no contribu- 
tion to the energy gap. We could have ignored it in Eq. 



(|26|l . In fact, the gap defined using Eq. i|22|) and taking 
e 2 /hi£i = e 2 / 'k£q is the so-called neutral energy gapi£ and 
it is equal to the other two gaps that we introduced in 
this section. 

Eq. ij3*§|) can also be written as 



A = 2E+ + [2e CS p + X (0)] 



e 

k£ q 



(40) 



In the lowest Landau level, the energy gap at v = 
1 is due to the excitation of a bimeron-antibimeron 
pair and the energy per electron in the UCS is 

e ucs {d) = [-*- \ [x(0)+X d (0)]] where X (0) = 

X R , L (N,N,0). Eq. igDJl can then be written, for this 
special case, as 



A = 2E+ + 2 [e ucs (d) - e UCS (d=0,t= 0)] 
which is just the form we used in Ref. 0. 



kIq) ' 
(41) 



V. INTERACTION BETWEEN 
QUASIPARTICLES 

With the simplifications introduced in the preceding 
section, the energy esc that enters Eq. (|22|l and Eq. 
is given by 



£SC 



•It 



Re[(p R > L (0))] 



(42) 



^EE^( G )(^' (-G)X^(G)) 

i,j G#0 



i,j G 

where, to simplify the notation, we have left implicit the 
index N of the Landau level. The soliton crystal is a su- 
perposition of a CSP with order parameters {(a* 1 - 7 (G)) } 
(computed at v — AN + 1) and a pure soliton crystal 
(PSC) with order parameters {(/3 lJ (G))} such that 

(p* (G)> = (G)) + (f* (G)> . (43) 
If we insert this decomposition into Eq. I|42|l . we find 

£sc — £csp fa + £csp-psc + £psc, (44) 

where 



2t 



scsp fa = -4Re[(a fl ' L (0))] 



(45) 



2 ^EE^( G ) (« M (-G)X^(G)) 

iEE^( G )(« 4J (- G ))(^( G )) 



i,3 G 
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is the energy per electron of the CSP (i.e. £csp i v ) 



2v 



£psc = -^Ee[(/3 R ' L (0))] (46) 

^EE^( G )(^(- G )>(^( G )> 

-4ee*« (°) (- g )> 

is the energy per electron of the PSC and 

£csp-psc = (47) 
+^EE H ^ («) (« Vi (-°)) (G)) 

1,3 G^O 

- = E E (G) (a^ (-G)) (G)) 

^ 1,3 G 

is the interaction energy (per electron) between the CSP 
and the PSC. 

The contribution epsc causes problem because it con- 
tains not only the energy to create the N qp solitons 
but also the interaction energy between the solitons. 
This interaction energy goes away in the limit Av — > 
0. As we said, however, we cannot go to arbitrarily 
small Ai/ numerically because solving the equation of 
motion for the single-particle Green's function then in- 
volves diagonalizing very large matrices. We must then 
find a way to remove the interaction energy in £psc- 
Two methods can be used. The first one is to re- 
place epsc by e PS c - £%nt where e int is the Madelung 
energy of the crystal of charged quasiparticles, assum- 
ing the quasiparticles to be point particles^. We re- 
fer to this method as the "Madelung" method. In the 
limit Av — * 0, the quasiparticles are very far apart 
and, if they have an isotropic charge distribution, it 
is a reasonable approximation. In the second method, 
which we refer to as the "form factor" method, we 
completely replace epsc ({(/3 iJ (G)/}) by the energy 
N qp e PSC (q))}) where e PSC ({(?£ (q)») is the 

energy per electron of a "crystal" of only one quasipar- 
ticle. In the case of solitons, which are quite extended 
and highly anisotropic objects it is necessary to use this 
second approach. 

To evaluate epsc {{{Pqp we ma ke use of the 

fact that, when the quasiparticles are very far apart (limit 
v — > 1) so that there is no overlap of the density or spin 
texture due to different quasiparticles, then we may think 
of the order parameters in real space as given by 



(/3«(r)> = ^^(r-R) : 



R 



where R is a lattice site. We know that 



(P* (r)> = 7F E (G)} ■ 

G 



(48) 



(49) 



but it is not possible to get hij (r) from this equation. 
We must make an approximation. Since we work in the 
low-density limit for the quasiparticles, it is a good ap- 
proximation to assume that for a "crystal" of one quasi- 
particle 

/ gi,j (r)) / V Eg (G)) e- Gr , if r e v c 
{P [r)) ip-\ 0, if v£v c ' 

(50) 

where v c is the volume of the unit cell centered at r = 0. 
Fourier transforming Eq. (|5U|I . we have 



(P iJ (q)> s 



N, 



dve^{(3^{v)) qp (51) 
L^(^3 (G) )A(q-G), 



G 



where the form factor 



A(q-G) = — { dre lqr e- ,;Gr . 

v c J Vc 



(52) 



depends on the shape of the unit cell of the soliton crys- 
tal. 

It now remains to compute the Hartree-Fock energy 
corresponding to the density and pseudospin textures 
given by the (/3jj (q)) ?p 's. The energy is still given by 
an equation similar to Eq. (|46[) where the summation 
W Eg i s now replaced by ^ Y^q- To go from the sum 
to the integral, we use 



2v ^ V ' 2v 



dq 



2v J (2tt) 2 
2nN e f dql 2 



(...) 



(53) 



2^ 2 J (2^) 2V '" ; ' 

Aso, because (0)) ~ 1/N V , we introduce a new 
field ,J (q) by the definition 



0^ (q) = N v (q)) 
With this last definition, we have 



QP 



(54) 



N, t e PSC ({(ft> (q)») (55) 

As a test of our "form factor" method, we have com- 
puted the energy gap due to the creation of bimeron- 
antibimeron pairs at v = 1 in the lowest Landau level 
N = 0. Figure |5] shows the energy gap computed from a 
triangular lattice of bimerons at v = 1.02 and t — 0.0025. 
In this case, the Madelung and form factor methods give 
identical results at small interlayer distances while the 
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Madelung method slighlty overestimates the energy gap 
at higher distances. The difference between the two ap- 
proches at large d is due to the fact that the charge 
density profile of the bimeron becomes more and more 
anisotropic as d increases. Also, the Coulomb interac- 
tion is stronger between point particles than between ex- 
tended particles so that the Madelung approach overes- 
timates the gap energy. 



1.2 



1.1 



Q_ 



co 0.9 

CT> 
CD 



CO 

O 



0.8 



0.7 



0.6 



0.5 



Supercell (Madelung) 
Supercell (form-factor) 
HF result 




0.2 0.4 0.6 0.8 1 

Interlayer distance (d/£) 

FIG. 5: The energy gap due to the excitation of a bimeron- 
antibimeron pair v = 1 computed using the form factor or 
the Madelung method and compared with the Hartree-Fock 
energy gap to the excitation of an electron-hole pair. 



To check the convergence of the supercell approach as 
the lattice constant gets very large, we show in Fig. 
the energy gap of the UCS at v — 1 computed at differ- 
ent values of v from a crystal of bimerons. The different 
curves in this figure are for different values of the tun- 
neling strength. The real gap of the system is, of course, 
defined for v — ► 1. We see that the gap converges more 
rapidly to its v — > 1 value when the tunneling is stronger. 
This is understandable since the size of a bimeron de- 
creases when t increases and, for sufficiently strong t, 
this size is independent of the lattice constant even at 
relatively high v. For smaller t the gap converges to its 
v — ► 1 value, but only at lower filling v. In the appli- 
cation of the supercell technique to the soliton gap in 
the next section, we will use the form factor method to 
remove the interaction energy. As we have just shown, 
this method is more appropriate in the case where the 
quasiparticle is highly anisotropic in shape. 




1.02 1.04 

Filling factor (v) 

FIG. 6: Energy gap of the UCS at v = 1 computed by the su- 
percell approach using the form factor method. The different 
curves are for different values of the tunneling strength t. 



VI. NUMERICAL RESULTS 

We now discuss our numerical results for the energy 
gap of the CSP. Our calculations are done in Landau 
level N = 2 around v = 9 using the form factor method. 
Figures|7Ja)-(c) contain our main results. Differents gaps 
are plotted as a function of the interlayer distance for 
tunnelings (a) t= 0.007; (b) t= 0.01; and (c) t = 0.02. 
The filled line is Ajjcs, the energy needed to create an 
ordinary electron-hole pair from the coherent liquid state 
at v = 9. At v = 9, the liquid state is unstable for d > d\ 
where the coherent striped state is the ground state. The 
Hartree-Fock gap represented by the curve with the filled 
squares is given by the energy to create an electron-hole 
pair in a coherent channel (see Fig. [21 where this gap 
is defined). The other curves give the energy gap calcu- 
lated in the supercell method for different filling factors 
v and the energy gap calculated with the field-theoretic 
approach explained in the appendix. 

From Fig. [7| it is clear that, in the CSP, the energy 
needed to create a soliton-antisoliton pair is smaller than 
that needed to create an electron-hole pair for typical 
experimental values of the tunneling parameter t. The 
transport gap is thus determined by the creation of these 
topological excitations (as it was the case for skyrmions in 
quantum Hall ferromagnet at v — 1 or with bimerons in 
bilayer quantum Hall systems)* 13 Figures a)- (c) show 
a rapid decrease of the energy gap near the transition 
between the coherent liquid and the CSP that should 
be observable experimentally. The curves corresponding 
to different filling factors show that the convergence of 
the supercell method is quite good near the liquid-CSP 
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0.8 



0.6 



0.4 — 
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t=0.007 



UCSHFGap 
CSP HF Gap 
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Supercellv = 9.0175 
Gradient approximation 



0.5 



1.5 



UCSHFGap 
CSP HF Gap 
Supercell v = 9.02 
Supercell v = 9.0225 
Supercell v = 9.025 
Gradient approximation 




UCSHFGap 
CSP HF Gap 
Supercell v=9.02 
Supercell v=9.0225 
Gradient approximation 




0.5 1 1.5 

Interlayer distance (d/t) 

FIG. 7: Different energy gaps in the UCS and CSP calculated 
as a function of the interlayer distance d/l and for different 
values of the tunneling parameter. For the supercell method, 
the gap is evaluated at different filling factors to show the 
convergence of the results to the true gap at v = 9. The 
gradient approximation refers to the field-theoretic method. 



transition but slow at larger values of interlayer distances. 
This slow convergence is due to the fact that the size of 
the soliton increases with interlayer distance as shown 
in Fig. [8] and the shape of the soliton is then restricted 
by the lattice constant as we explained previously. As 
d/£ increases, it becomes necessary to go to lower filling 
factors to achieve convergence, something we cannot do 
numerically. In any case, the soliton gap is always lower 
than the Hartree-Fock gap at higher values of d/£ since 
our approach overestimates the energy gap. Increasing 
t decreases the size of the solitons, however, so that it 
is possible to achieve better convergence by increasing 
the value of the tunneling parameter t. This is seen by 
comparing Fig. [3(a), (b) and (c). Notice also that, for 
smaller solitons, the soliton gap is closer to the Hartree- 
Fock result, as expected. 




1 1.2 1.4 

Interlayer distance (d/t) 

FIG. 8: Soliton size calculated with the supercell (filled sym- 
bols) and field-theoretic (empty symbols) methods as a func- 
tion of the interlayer distance at v = 9.02. In the supercell 
approach the size of the soliton is found by fitting the y depen- 
dence of the phase in a channel with ip(y) =4 tan" 1 |^e~ !/ ' /Ls j . 

We also show in Fig. [7| the gap calculated with the 
field-theoretic method (see Eq. (|17fl b This gap has the 
same qualitative behavior with interlayer distance, ex- 
cept at small d near the phase transition. It is larger 
than the gap calculated in the microscopic approach. As 
we explain in the appendix, the field-theoretic result is 
incorrect at small d or large t (fig. Wi c )) where the stripes 
are not fully developped. At large ri, we cannot say how 
different the two gaps (macroscopic and field-theoretic) 
are because the gap found in the microscopic approach 
has not yet converged at the lowest filling factor we can 
achieve. 

In the field-theoretic method, the soliton size, L*, is 
obtained by the procedure outlined in Sec. Ill . When 
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the Coulomb interaction between parts of the soliton is 
properly included, we find numerically that L* s increases 
with d as in the supercell calculation. Both approaches 
give the same trend for the soliton length. The detailed 
behaviour with d/£ is quite different, however. Ccarly, 
the field-theoretic calculation does not capture all the 
subtleties of the We recall that, as the interlayer distance 
increases, the width of the LCR's becomes smaller. The 
behavior of the soliton size may be understood as aris- 
ing from the Coulomb energy, which favors spreading the 
charge of the soliton. Our results are plotted in Fig. [S] 
In this figure, we see that the supercell and field-theoretic 
results do not match for large t. This is again due to the 
fact that the stripes are not fully formed at large t so that 
the expression of Eq. (|A.32J| for the topological charge is 
not correct. As expected, Fig. [S] shows that the soliton 
size decreases with t. 

We have neglected quantum fluctuations in our calcu- 
lation. These fluctuations increases in importance as d/i 
increases. They renormalize the pseudospin stiffness and 
will probably also change the size of the solitons and the 
quantitative values of the energy gaps. Inclusion of these 
fluctuations is, however, beyond the scope of this paper. 

VII. CONCLUSION 

We have computed the energy gap due to the creation 
of a soliton-antisoliton pair in the linearly coherent chan- 
nel of the coherent striped phase found in higher Landau 
levels in a bilayer quantum Hall system. We have com- 
puted this gap using a microscopic unrestricted Hartree- 
Fock approach as well as a field-theoretic approach valid 
in the limit of slowly varying pseudospin texture. With 
both methods, we find that the this energy gap is lower 
in energy than the Hartree-Fock gap due to the creation 
of an electron-hole pair in a coherent channel (a single 
spin flip) so that solitonic excitations must play an im- 
portant role in the transport properties of the coherent 
striped phase. 
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APPENDIX: MICROSCOPIC EXPRESSIONS 
FOR THE PARAMETERS OF THE 
FIELD-THEORETIC MODEL 

In this appendix we present the details of the deriva- 
tion of the microscopic expressions for the parameters p s 
and T used in the field-theoretic model of Sec. III. We 



drop the Landau level index N here since all order pa- 
rameters are to evaluated in the partially filled level N. 
We begin by defining the pseudospin density operators 

p(q)=P^(q)+/' L (q), (A.l) 

PM = 1[P R ' R (<1)-P L ' L (<1)], (A.2) 

^(q) = ^[p R ' i (q)+p I " fi (q)], (A.3) 

P,(q) = ^[p fl ' i (q)-p L ' fl (q)]- (A.4) 

The total Hartree-Fock energy of the electrons in the par- 
tially filled level for an unbiased bilayer can be written 
as 

E HF = e (£) , (A.5) 

where 

e = -2N^t( Px (0)) (A.6) 

+ljV^T(q) (p(- q ))( p(q) ) 
q 

+A^J 2 (q) ( Pz (-q))<p, (q)) 
q 

+iV^Jx(q) [( Px (-q))(p, (q)) 
q 

+ <M-q)> (py (q))]' 

We have introduced the interactions 

J z (q) = H RtR (q) - H RtL (q) - X R>R (q) , (A.7) 



T (q) = H R , R (q) + H R . L (q) - X R , R (q) , (A.8) 

and 

J X (q) = -X R>L (q) . (A.9) 

In Eq. (|A~6jl . H R , R (0) = H R:L (0) = because of the in- 
teraction between the 2DEG and the positive background 
of the donors. 

We now introduce a unitless and unitary pseudospin 
field S a (r), with a = x, y,z related to the guiding cen- 
ter density operators in the pseudospin formalism by the 
relation 

S a (r)=4rf 2 A^(p Q (r)), (A. 10) 

and a projected^ electron density by the relation 

n(r)=A^(p(r)). (A.ll) 
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Using the definition of the pseudospin operators S a (r) 
and taking the Fourier transform of Eq. (|A.6|I . we have 



a channel is given by 



-t 

1 

+ 8^ 2 



dvS x {v) 



(A.12) 



l^Jdv f dr'j ± (r-r')S ± (r)-S±(r') 
^zjdrj dv'S z {v)J z {v-v')S z (v') 
J dv J dr'ra(r)T(r-r>(r'). 



Writing S a {r) in spherical coordinates, it is easy to 
describe the CSP ground state as 



S x (r)csp = sm6(x), 
Sy(r) CS p = 0, 
S z {r)csp = cos9(x), 



(A.13) 
(A.14) 
(A.15) 



while the density (n(r)) = est is uniform. For a state 
where there is a spin texture only in the channel centered 
at x = (channel 0) while the other channels remain in 
their CSP ground state configuration (we recall that £ is 
the interstripe distance), we write 



S x (v) 



sin 9(x) cos (p(y), if |x| < |, 



sin 8{x), 



if N > f, 



(A.16) 



S y (r) 



sin 8(x) sin <p(y), if |i| < |, 
0, ifN>|, 



S z (r) = S z (r) CSPl 



n(r) = n(r)csp + Sn(r). 



(A.17) 



(A.18) 



(A.19) 



In these equations, 0(x) is given by its value in the CSP. 
Defining 



Ji,j{y~y')= dx dx'J±(r — r') sm8(x) sm0(x'), 

(A.20) 

where Ci corresponds to the i-th channel of width £/2 
centered at Xi and J c — f x l ^u > ^ is easy to show that 
the energy difference between the this last state and the 
CSP ground state i.e. the energy to create one soliton in 



fe = 



-t 
2^P 



dxsm8(x) / dy [cos ip(y) — 1] 



4,^2^1 / dy J ' d y' J *,o{y-y') [cos^(y') - 1] 



dy / dy'Jofiiy - y') [cos(y>(y) - <p{y')) - 1] 



4-/*/ dv'5n{r)T{r - r')fc(r') 
J dv J dr'Sn(r)T(r -r')n(r')csp- 



(A.21) 



The first two terms contribute to the effective tunnelling 
term T while the third term is directly related to the 
pseudospin stiffness of the system. The fourth term takes 
into account the Coulomb interaction between different 
parts of the soliton and the last term is the interaction 
between the charge of the soliton and that of the CSP. In 
an antisoliton, this fifth contribution would have exactly 
the same value but with opposite sign so that this last 
term does not contribute to the transport gap. 



1. Calculation of the pseudospin stiffness p s 

To extract the pseudospin stiffness from the third term 
of Eq. (|A.21|I , we make a long- wavelength expansion of 
the cos(ip(y) — ip(y')) — 1 term. This expansion is possible 
if the pseudospin texture varies slowly in comparison with 

Jo,o(y)- We S et 



dy / dy'j 0t0 (y - y) [cos(y(y) - <p{y')) - 1] 



8nP 



1 



dy' y' 2 Jo,o(y') 

j v u- y j 

(A.22) 

Comparing this last result with Eq. (112H , we see that 
1 



8ir£ 2 



dy y Jo,o{y)- 



(A.23) 



The pseudospin stiffness can be written, more explicitely 

as 

^ = -8^/^I^£<MqK- (A.24) 

J q 

x j dx j dx' sm8(x)sme{x')e lq * {x - x '\ 

J Co J Co 

with L x and L y the length and width of the sample. This 
allows the integrals over x and x' to be totally decoupled. 
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In fact, defining the form factor 



c 



dxs\rL9(x)e lq * 



(A.25) 



£^(p*(Gx)) 



sin[(G x - g:c )g/4] 
(G x -q x )£/4 



we can write 
1 



Ps 



dq x \Vt{q x 



|2 <* 2 J ± (q) 



(A.26) 



The form factor O^) takes into account the influence 
of the shape of the charge modulation along the x axis 
in the CSP phase on the effective pseudospin stiffness in 
the one dimensional sine- Gordon model. 



2. Calculation of the tunneling parameter T 



sum over the reciprocal lattice vectors of a ID lattice of 
lattice constant £/2, noted G x , and 

J d v J d y'Ao(y~y')[cosip{y')~i] 

= i^J ± (G ;c ,0)|fi(G :c )| 2 -i-^^J i (g :c ,0)|^)| 2 . 



Combining the two terms, we find 



(A.29) 



T 



2ttP 



~^^2J±(Q x ,o)\n(q x 



The effective tunnel coupling T can be extracted from 



the first two terms of Eq. (|A.21|> . The first term renor- 
malizes the tunnel coupling in the ID effective theory, 
taking into account that interlayer coherence exists only 
in the LCR's. This first term is simply 



2ir£ 2 



O(0) / dy [cos <p(y)- I}. 



The second contribution to the effective tunnel cou- 
pling comes from the exchange energy between channel 
(where a pseudospin texture was created) and the other 
channels. In these other channels, the in-plane pseu- 
dospin component is totally polarized along the x di- 
rection and the exchange interaction between channel i 
and channel favors a configuration in channel where 
the pseudospin is also polarized along +x, just like the 
simple tunnel coupling t. In other words, there is an en- 
ergy cost, even in the absence of tunneling, to make a 
rotation of the pseudospins in one channel because of the 
interaction with the pseudospins in the other channels. 

It is possible to extract a simple form for this coupling 
from the second term of Eq. I|A.21(1 



since 



2 e iq x { Xi - Xo ) 



i#0 



—no) 



£ J dyJiAv) = ^^Ji(fe,0)|% 

x i#0 gx 

= ^££Jx(&,0)|fi(<fe)|V*-<* 

X i <?x 

-^-^J±(fe,0)|O(g x )| 2 , (A.28) 



with x n — Xq = n£/2 the center-to-center distance be- 
tween channels n and 0. Because there is a sum over the 
channels, the sum on the wave-vectors q x reduces to a 



3. Sine-Gordon soliton and the Coulomb energy 

If we combine the tunneling and exchange terms, we 
find that the energy cost to make one soliton localized 
in a channel of the CSP is given by Eq. 112(1. As we 
mentionned in Sec. Ill, the static solution that minimizes 
this energy functional is the sine-Gordon (or kink) soliton 



(A.27) tp(y) = 4t&n 



We now add to Eq. i|12|) the Coulomb interaction en- 
ergy between different parts of the soliton 



SE, 



Coul 



J dv J dr'Sn (r) T (r - r') Sn (r') . 



(A.31) 

To relate Sn(r') to the angles 8 and ip, we use the def- 
inition of the topological charge density given in Eq. 
pi[l. We assume that, in the one-soliton state, only ip (y) 
changes along a channel and that 9 (r) is given by its 
value in the CSP. We have 



Sn (r) 



- — Vp(r) x (Vcos0(r)) • z 

47T 



(A.32) 



1 d<p(y) d 
Ait dy dx 



cos 9 (x) 



At this point, we must remark that if we use 
the sine-Gordon solution in Eq. (|A.32|) and integrate 
the projected density Sn (r) in a channel, we find 
J-£/4 dx J-m dySn (r) = 1 only if cos 9 (x) varies from 
— 1 to +1 in the channel i.e. only in the limit or large in- 
terlayer distances where the stripes are fully developped. 
In consequence, we do not expect our field-theoretic 
model to be valid near the transition between the UCS 
and the CSP. 

We insert Eq. l|A.32p into Eq. HA.31J1 . an d define the 
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form factor (for a channel centered at x = 0) 

A(q x ) = [ dxe- iq * x -^cos6(x) (A.33) 

JCn dx 



sin(q x ~ G x )£/4 



(fe-Gx)C/4 ' 
and the effective interaction V e g (y — y') in a channel 

Veff (y-y') = W \ A (^)l 2 T W eM*"*'). (A.34) 



We then find for the Coulomb interaction 
oEcoui = T^zt I dy dy — — V e g (y - y ) 



32tt 2 



dy 



dy' 
(A.35) 

If we add the contribution SEcoui to Eq. 1121) and min- 
imize the energy with respect to (p (y) , we find that it 
introduces a nonlocal term to the sine-Gordon equation. 
The resulting equation is then very difficult to solve. To 



get an approximation for the Coulomb energy, we de- 
cided to proceed in the following way. We take, as a trial 
solution, the kink soliton 



<p(y) = 4 tan 



,-y/K 



(A.36) 



where L* is the width of the soliton. The Coulomb energy 
is then 



SEcoui (L*) 



tt£ 2 
32^2 



dq\A(q x )\ 2 T(q) sech : 



The total energy for the soliton is 

Ps 



E 



4TL* + 5E 



Coul 



( Kq v L* s 

V 2 
(A.37) 



(A.38) 



We find L* by minimizing numerically the total energy 
E. In our numerical calculation, we use T (q) = Hm (q) 
instead of Eq. (|A.8(I . This is also the interaction consid- 
ered in similar calculations^*^. The use of Eq. 
leads to non-physical results. 
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